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Abstract 

m 
o . 

We study the radial pulsation frequencies of slowly rotating neutron stars in general 



relativistic formalism using realistic equations of state. It is found that the pulsation 
frequencies are always an increasing function of rotation rate. The increasing rate of 



CZ3 

frequency depends on the nature of equations of state. 



Cameron suggested that the vibration of neutron stars might excite motions that might have 
interesting astrophysical applications, which lead to a series of investigations of the vibrational 
properties of neutron stars. The earliest detailed calculations were done by Meltzer and Thorne 
0] and Thorne ]3|], where they investigated the radial as well as nonradial oscillations using 
available equation of state, such as the Harrison - Walker - Wheeler equation of state. These 
and other early studies by Wheeler Q], Chau f| and Occhionero |J indicated that the majority 





of the fundamental mode radial oscillation periods for the neutron stars would lie in the vicinity 
of about 0.4 msec, and, that the first few quadrupole oscillations would have periods that are 
also fractions of a millisecond. Moreover, these oscillation periods were estimated to be damped 
by gravitational radiation with damping timescales of the order of one second. 

If an oscillating neutron star is also rotating, then rotation can provide a coupling between 
the radial and the nonradial modes, which leads to rapid loss of radial vibrational energy H. 
As pulsars are supposed to be rotating neutron stars, the study of radial pulsation frequency 
of rotating neutron stars constitutes an interesting astrophysical problem to investigate. The 
problem in its details has not been worked out yet in a fully general relativistic framework 
with arbitrary rate of rotation, as it is a very complicated task. In recent years, a Newtonian 
theory for the above purpose became available mainly due to the work of Ipser and Managan 
PJ Managan l| and Ipser and Lindblom 11. More recently, a formalism for computing the 



oscillation frequencies of rapidly rotating neutron stars is given by Cutler and Lindblom (1C 
using the post-Newtonian approximation. They employed polytropic equations of state, which 
are an idealization to describe neutron star matter at high densities and which are unlikely to 
reflect accurately the adiabatic index as function of density, an important ingredient to describe 
the stability of oscillation frequencies. Therefore, it is desirable to use more realistic equations 
of state for neutron star matter and a general relativistic formalism. 



Chandrasekhar and Friedman |TT] have given a general relativistic formalism to calculate 
the effect of rotation ( to order Q 2 , where Q is the angular velocity of rotation ) on the eigen- 
frequencies of radial pulsations of stars. This formalism is an exact formula to calculate the 
frequency (cr') of oscillations of a 'slowly' rotating stellar configuration which depends only on a 
knowledge of the Lagrangian displacement associated with the radial mode of oscillation of the 



nonrotating configurations and of the uniform spherical deformation caused by the rotation. 
In this paper, we apply this formalism using recently suggested realistic equations of state to 
estimate the rotationally shifted eigenfrequencies of radial oscillations of neutron stars. 
The Chandrasekhar-Friedman formula is of the form: 

a' 2 h = I 2 + I 3 + h, (1) 

where 
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The subscript (S) refers to Schwarzschild (i.e., nonrotating configuration) values and w(r) is 
the angular velocity of the fluid element relative to the local inertial frame. 

In the above equations only the spherical deformation effects of rotation are manifest. Con- 
tributions due to the quadrupole deformation are all zero. The metric that will correspond to 
this was given by Hartle [|12j (see also Hartle and Thorne [13| 



ds 2 = -e u c 2 dt 2 + r 2 sin 2 e(d(t) - wdt) 2 + e x dr 2 + r 2 d6 2 (6) 
where u(r) is the angular velocity of the cumulative inertial frame dragging. In terms of uj(t), 



the angular velocity of the local inertial frame is written as 

Co = Q — u, (7) 

where Q is the angular velocity as seen by a distant observer. In Eq.(|6|), one can write 

v = V(s) + U(n), (8) 
A = A {S ) + A(n), (9) 

to distinguish the corrections of 0(Q 2 ) due to the rotation from the Schwarzschild values by 
the subscripts (O) and (S) respectively. In Hartle's notation, 

= h (10) 

and 

\ n = m ° (11) 

r — 2Gm(r) / c 2 ' 

where h and m Q are functions of r. To 0(Q 2 ), 

v /=^ = r 2 sin6[l + (A + v) { n)}e {x+u) ^ , (12) 

u° = (l- m + ±V 2 )e-^, (13) 

V = ru)sin9e~^ s \ (14) 

and 

u\ = rVsinO. (15) 
The various other quantities appearing in Eqs. (@ - 0) £1X6 cLS follows: 
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The equation governing infinitesimal radial pulsations of a nonrotating star in general rela- 
tivity was given by Chandrasekhar WM, and it has the following form : 



where £(r) is the Lagrangian fluid displacement and ca is the characteristic eigenfrequency (c is 
the velocity of light). The quantities F, G, H depend on the equilibrium profiles of the pressure 
and density of the star and are given by 

F = -e- x e u Tp/(p + pc 2 ), (21) 

„ A ,. r_ /l dv ldX 2\ 

G = -e x e u \Tp( — - H + - + 

I- y \2dr 2dr r> 

p? + r ?l/(p + pc 2 ), (22) 



dr dr ■ 

p + pc 2 l r dr p + pc 2 ' c 4 
where T is the adiabatic index, defined in the general relativistic case as 

r = ( 1 + <^5(^)' (24) 



and 
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Note that in Eqs.( pT| - |25|) , p(r), p(r), v(r) and A(r) correspond to the values in the nonrotating 
configuration, whose profiles are obtained from the equations for the hydrostatic equilibrium of 
nonrotating degenerate stars in general relativity |L5 



dp G(p + p/c 2 )(m + Anr^p/c 2 ) 
dr r 2 (l — 2Gm/rc 2 ) 



(26) 



dm 9 , _, 

- = W„ (27) 



du 2G (m + 4vrr 3 P/c 2 



(2f 



dr r 2 c 2 (1 — 2Gm/rc 2 ) ' 

A = - ln(l - 2Gm/rc 2 ). (29) 

The boundary conditions to solve the pulsation equation (|20"D are 

C(r = 0) = 0, (30) 

dn e v l 2 8 

5p(r = R) = - Tp— Y ye-l%\ r = R = (31) 

Since p vanishes at r = R, it is generally sufficient to demand 

£ finite at r = R. (32) 



Equation (p0[), subject to the boundary conditions Eq. (p0|) and (pi]), is a Sturm - Liouville 
eigenvalue equation for a 2 . The following results follow from the theory of such equations : 

1. The eigenvalues a 2 are all real. 

2. The eigenvalues from an infinite discrete sequence 

a 2 < of < . . . < a 2 n < , (33) 

with the corresponding eigenfunctions £o( r ) ; £i( r ); ••• ) £n( r )- The importance of this 
result is that if the fundamental radial mode of a star is stable [o 2 > 0) , then all the 



radial modes are stable. Conversely, if the star is radially unstable, the fastest growing 
instability will be via the fundamental mode (a 2 more negative than all other cr^). 

3. The eigenfunction £ D corresponding to a 2 Q has no nodes (fundamental mode) in the interval 
< r < R; more generally, £ n has n nodes in this interval. 

4. The £ n form a complete set for the expansion of any function satisfying the boundary 
conditions Eqs. fl3"0|) and (PH). 

Given an equation of state p(p), Eqs. (p6|) and (|29|), called Tolman-Oppenheimer-Volkoff (TOV) 
equation, can be numerically integrated for a given central density to obtain the radius R and 
gravitational mass M = m(R) of the star. 

For the rotational spacetime corresponding to the metric Eq.(|), the stellar structure equa- 
tions fl26|) - fl29f) must be supplemented by additional equations that provide the profiles for 
ijj(r) and the (rotation induced) pressure and mass deformation terms. The details of which are 



given in references [[T2|] and |[L3|| , and so are not reproduced here. For the integrand appearing 
in Eq.(|4]), the subscript (S) refers to the Schwarzschild values of quantities p, p etc. In the 
Eqs.(0), ([jD and (|5|) all the quantities must incorporate the rotational contributions. 

The main ingredient to solve for the structure and pulsation features is the equation of state 
(EOS), p = p(p). The range of p expected for neutron stars spans a very wide range: from 
about 7.8 g cm~ 3 near the surface to p > 10p o in the interior, where p Q is the equilibrium nuclear 
matter density (=2.8 xlO 14 g cm~ 3 ). Structure parameters of a neutron star depend sensitively 
on the EOS at high densities [IB, |T7| especially around the density region 10p o . Therefore, the 
oscillation features are also expected to possess marked sensitivity on the high density EOS, 
and so it becomes imperative to use realistic neutron star EOS model rather than a polytropic 



form for the EOS. 

Despite two decades of theoretical research, there exists no clear consensus on the behaviour 
of the EOS at high densities. The nearly two dozen EOS models available in the literature 
comprise a rather broad set. For our purpose here, we adopt a realistic non-relativistic model 



given by Wiringa, Fiks and Fabrocini |L8| and a relativistic model given by Sahu, Basu and 



Datta [19|. The non-relativistic model describes the high density matter in the interior region 
of neutron star as beta-stable neutron matter. This model possesses desirable properties for 
two-nucleon scattering phase shifts, equilibrium nuclear matter and binding energies of light 



nuclei. The three versions are given by Wiringa, Fiks and Fabrocini |18[ which correspond to 
a variation in the choice of the nuclear interaction involved. We choose here the UV14+UVII 
(neutron, proton, electrons and muons). Whereas, in recent times, the relativistic model has 
drawn considerable attention. In the relativistic approach, one usually starts from a local, 
renormalizable field theory with baryon and explicit meson degrees of freedom. The theory is 
chosen to be renormalizable in order to fix the coupling constants and the mass parameters by 
empirical properties of nuclear matter at saturation. As a starting point, one chooses the mean 
field approximation which should be reasonably good at very high densities. This approach 
is currently used as a reasonable way of parameterizing the EOS. However, in recent years, 
the importance of the three-body forces in the EOS at high densities has been emphasized 



by several authors (Jackson, Rho and Krotscheck f20 |; Ainsworth et al. |2T|). This gives 
theoretical impetus to study the chiral sigma model, because the non-linear terms in the chiral 
sigma Lagrangian can give rise to the three-body forces. In the present calculation, we take 
the SU(2) chiral sigma model EOS for the neutron-rich matter in beta equilibrium fl9f . The 
composite equations of state to determine the neutron star structure was set up as follows. For 



the region 3.7 x 10 11 < p < I0 14 g cm 3 , p(p) is taken from Negele and Vautherin [22|, while 



for p < 3.7 x 10 g cm , p(p) is taken as given by Baym, Pethick and Sutherland |23| and 



Feynman, Metropolis and Teller \24 \. 



We first solved the stellar structure equations that provide the nonrotating neutron star 
configuration. The numerical integration was done using the modified Euler method with 
variable step size (where density gradient is large) and a logarithimc interpolation for the 
(p, p) data. The method was tested by requiring agreement with the results for the neutron 



star mass, radius and momemt of inertia as reported by Wiringa, Fiks and Fabrocini JT8 
The rotational contribution to p, p, m and v were calculated similarly (to accuracy < 0.6%, 
compared with the published results by Hartle and Thorne |0| for the case of Harrison- Walker- 
Wheeler EOS). The nonrotating general relativistic pulsation Eq.(|C]) was replaced by a set of 
difference equations which were cast in tridiagonal form. The eigenvalue and eigenfunction 
were evaluated using EISPACK routine. The eigenvalues calculated by us agree with the radial 



pulsation frequencies of nonrotating neutron stars for the Wiringa, Fiks and Fabrocini [^] EOS 
calculated by Cutler, Lindblom abd Splinter Finally, the integrals Ji, I2, h and 1^ were 
evaluated using standard integration package, and then substituted in Eq.([TJ) to obtain the 
value of a', the eigenfrequency of raidial pulsation of the star in rotation. This was done for 
several values of the rotation rate (f2) to obtain a'(fl). For purpose of illustration, we consider 
only one configuration, namely a 1.4 M & (M = solar mass) neutron star. Masses of neutron 
stars estimated from binary pulsar data are all consistent with a value 1.4 M e . 

The results of our calculation are presented in figures 1 and 2 and table 1. In figure 1, we 
show the range in the values of the frequency of radial oscillation for the nonrotating neutron 
stars configurations. In this figure, P stands for the period of oscillation and is plotted as 
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a function of the stable gravitational mass of the nonrotating neutron stars (with maximum 
stable mass equal to 2.59 M for EOS |D] in Fig. 1(a) and 2.18 M© for EOS 0] in Fig. 



1(b)); curve 1 refers to the fundamental mode of oscillation and the curves 2 - 5 to the first four 
harmonics. For the fundamental mode, we find that P lies in a rather narrow range of values: 
(0.3-0.5)msec for both the EOS. For the higher modes considered in Fig. 1(a) and 1(b), values 
of P lie roughly between the limits 0.25 msec and 0.1 msec, and seem to saturate at about 0.1 
msec. In case of quark stars [56], the fundamental mode of oscillation periods are in a range of 



values (0.05 - 0.3) msec and for higher modes, the periods are < 0.1 msec. Thus, the neutron 
star periods are slightly higher than for quark stars. The pulsation frequency as a function 
of fl/Q c for the fundamental model (n=0) for both relativistic EOS [|19|] and non-relativistic 



EOS [18 1 are shown in Fig. 2(a) and 2(b) respectively. In both the cases, it is found that the 
pulsation frequency is an increasing function of Q. In Fig. 2(b), rate of increase is higher, 
which is possibly due to the nature of interaction used in EOS. The EOS is stiff for relativistic 



case 



[Bj] , whereas, it is soft in the non-relativistic case [|Tj|. However, we do not find such a 
feature in the nonrotating case. Table 1 presents the values of a'(Q) for the fundamental mode 
and the first two harmonics for both relativistic EOS [TJSj (a) and non-relativistic EOS |TB[ (b). 
For comparison, we have included in this table the values of a, the pulsation frequency of the 
nonrotating stars. The parameters appearing in this table are all for a 1.4 M & neutron star. 
We have chosen the rotation rate Q to vary from 0.m c to 0.5fi c , where fi c = (GM/R 3 ) 1/2 . The 
upper limit 0.5Q C is almost equal to the termination point 0.6f2 o , where f2 = (3GM / '4R 3 ) 1 / 2 , 
for sequences of rotating stellar models, applicable to both Newtonian and general relativistic 
stellar models and is independent of the equation of state of the stellar matter [pTfl . We have, 



however, not investigated the dependence of this termination limit on the pulsation degrees of 

9 



freedom. 

From table 1, the following behaviour for a'(Q) emerges. For the fundamental mode (n=0), 
a 1 possess a monotonic dependence on Q and it is found to be always an increasing function of 
Q. The same feature is also found for n=l and n=2 modes for both EOS (a) and (b). In EOS 

(a) , the a' is very sensitive to Q for n=0 and n=l. This is due to the nature of force involved 
in EOS. On the otherhand, a' seems to saturate for n=2 in EOS (a) and n=0 and n=l in EOS 

(b) . 

We have thus applied the general relativistic formalism by Chandrasekhar and Friedman 



I I ] to calculate frequencies of radial pulsations of rotating neutron stars, for realistic EOS. 



However, the formalism is based on the use of Hartle's prescription |T_2| to treat the rotation 
general relativistically (to order f2 2 /f2 c 2 ). The essential physics of high density matter through 
the incorporation of a realistic EOS is important, and this feature is present in our work. The 
main conclusion of our work is that the magnitude of the effect of the dynamical coupling of 
oscillation and rotation will differ substantially in a general relativistic treatment of the problem 
as compared to Newtonian and post- Newtonian approximation approaches. Extension of this 
work to include the damping time scale with other different EOS will be reported in a future 
work. 

It is great pleasure to thank B. Datta and S. S. Hasan from Indian Institute of Astrophysics, 
Bangalore, for useful and fruitful discussions during the preliminary stage of this work. 
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Figure 1: Periods (P) of radial oscillations versus stable gravitational mass of nonrotating 
neutron stars: (a) relativistic EOS [19] and non-relativistic EOS [18]. Curve 1 is for the 
fundamental mode and curves 2-5 correspond to the first four harmonics. 



Figure 2: The pulsation frequency (a') as a function of Q/Q c of rotating neutron stars: (a) 
relativistic EOS [19] and non-relativistic EOS [18]. 
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Table 1: Pulsation frequency as a function of the rotation rate, <r'(Q), for n=0, 1 and 2 radial 
modes, a stands for the radial pulsation frequency of the nonrotating neutron stars. The 
neutron stars configuration correspond to a mass of 1.4M Q with: (a) relativistic EOS [19] and 
(b) non-relativistic EOS [18]. 



n 




(a) 




(b) 








a 


a'(n/n c ) 


a 


a'(n/n c ) 






(ioS- 1 ) 


(io^- 1 ) 


(ioS- 1 ) 


(lO 4 ^ 1 ) 





0.1 


2.108 


5.614 


2.066 


2.081 




0.2 




6.625 




2.123 




0.3 




8.654 




2.184 




0.4 




10.464 




2.251 




0.5 




11.839 




2.313 


1 


0.1 


3.666 


6.598 


4.979 


5.012 




0.2 




7.529 




5.102 




0.3 




9.062 




5.222 




0.4 




10.536 




5.340 




0.5 




11.396 




5.436 


2 


0.1 


5.379 


5.721 


7.155 


7.400 




0.2 




6.648 




7.963 




0.3 




6.540 




8.545 




0.4 




6.776 




8.975 




0.5 




6.847 




9.213 
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